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ABSTRACT 

An isothermal elastohydrodynamic lubrication of 
materials of both low and high elastic moduli was being 
Studied. A physical elastohydrodynamic lubrication model for 
materials of low elastic modulus was postulated. The 
governing equations for the model are the momentum and 
continuity equations for fluid and the Hertzian equations 
for solid. A numerical solution of finding the minimum film 
thickness of the form hj,;,/(U*)'’* was presented over a wide 
range of loads, between w*/(U*)'’? = 0.9 and 20. The results 
were found to have a good agreement with the earlier 
experimental results. 

A theoretical approach was also used to analyze the 
elastohydrodynamic lubrication problem under laminar regime. 
The governing equations for the theoretical analysis are the 
basic elasticity equation together with the momentum and 
continuity equations. Results of film thickness and pressure 
profiles for different hj;,, were obtained, which were’ then 
used tortand velocitywdarstribution, load capacity, and drag. 
The impacts of pressure-dependent viscosity and elastic 
deformation on the analysis of elastohydrodynamic 
lubrication were examined. The variations of load and drag 
with hg;, were noted. Present theoretical results were found 
to have a good agreement with the earlier experimental and 


theoretical results. 
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1. INTRODUCTION 


1.1 Statement of the Elastohydrodynamic problem 

Elastohydrodynamic lubracation, unlike classical 
hydrodynamic lubrication, takes into consideration the 
variation in the viscosity of a lubricant film with pressure 
and the elastic deformation of the bounding solids. The 
solution of film shape and pressure profile for a contact 
thus requires a Simultaneous solution of the continuity and 
momentum equations for the fluid and the elasticity equation 
for the solid with a viscosity which changes along the film 
as the pressure varies. 

When the bounding solids are made of materials of high 
elastic modulus es ga metal, elastic deformation is 
relatively small as compared to the film thickness. High 
film pressures are generated and cause a large changes in 
lubricant* viscosity along the film. On the other hand for 
low elastic modulus bounding solids ,e.g. rubber, elastic 
deformation is relatively large as compared to the film 
thickness. The film pressures generated between soft 
contacts are insufficient to cause significant increases in 
the lubricant viscosity ,in Sharp contrast with the metallic 
@aSeneuTheretorejevthe: blubracant in this: case may be 
considered isoviscous. 

Though,@much attenmitionBhas® beenspaid: uto\ cthe@pproblems 
concerning the hubpocarren eobeemetallbicm gearsysiroller 


bearings and cams, recently there is an increasing interest 
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in the lubrication of . soft elastic surfaces which finds 
veselfviimpontantuginesathevalubrication&® off  foritminstance, 
reciprocating seals, flexible-pad thrust bearings and motor 
car tyres aquaplaning on wet roads. 

In the present study, analysis will be made under 
laminar regime for both high and low elastic moduli bounding 
solids under isothermal elastohydrodynamic lubrication 
condition? 

The geometry of the bounding solids are chosen as two 
Rotating elastroscyhinderscwith a) uthin i lbbubricant vefilmh in 
between. This particular geometry is chosen because of the 
fact that many contacts between machine components can _ be 
represented by cylinders which provide good geometrical 
agreement with the profile of the undeformed solids in the 
immediate vicinity of the contact. 

In summary, results obtained from the present study are 
(i) £Plmeprnofiles. 

(ii) pressure profiles from which load capacity can be 
obtained. 
(ei h)mavé boc ipy eiprofhiles bef rempewhiche dragthioncesteane be 


computed. 


1.2 Review of relevant literature 

The early interest in the subject of elastohydrodynamic 
lubrication arose from the study of gear lubrication. Martin 
[1] approached the gear lubrication problem by considering 


the equivalent problem of the lubrication of a cylinder near 
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a plane with the assumptions that the lubricating contacts 
were rigid and lubricant was isoviscous and incompressible. 


From his work, three useful relations were obtained. 
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However, film thickness predicted by his theory could not 
account for the existence of fluid films which were thick 
enough to prevent severe metallic contact in gears. 

In the early 1960's, gears lubricated by a hydrodynamic 
action was confirmed by direct film thickness measurements 
in disc machine designed to simulate gear tooth contact. 
Crook [2], using a capacitance technique, measured the film 
thickness which was of the order of 2x10°4mm (10°° inches). 

Peppler [3] and Meldahl [4] examined tthe effect of 
elastic deformation in the gear lubrication problem. Neither 
Peppler's nor Meldahl's work Suggested that elastic 
deformation alone could be responsible for the existence of 
a continuous fluid film between gear teeth. 

Gatcombe [5] examined the agence of high pressure 
upon! lubricant: viscosity.*His work suggestedothatmthe effect 
Olmmtheieviceosity-pressure »characteristics® of sa, lubricant 
alone could not account for the existence of a continuous 


fluid film in gear contacts. 
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Grubin [6] developed an approximate film thickness 
formular which included the effect of elastic distortion and 
pressure dependent viscosity. He assumed that the elastic 
distortion would take the form of Hertzian deformation under 
highly loaded conditions. His predicted values of the film 
thickness were orders of magnitude greater than that of the 
Martin's, and consistent with the formation of a continuous 
fluid film between gear teeth. One useful relation emanating 
from his work in terms of dimensionless quantities is 
Has bis to Ga) Sob sew ibe t' (he2) 
However, the exact features of the contact mechanism do not 
emerge from thie Simple analysis. 

Dowson and Whitaker [7] presented a solution to the 
problem of elastohydrodynamic lubrication of highly loaded 
elastic cylinders with the variation of viscosity with 
pressure under isothermal condition. A numerical method is 
developed which enables a pressure curve to be found which 
Satisfies the elastic and hydrodynamic compatibility of the 
System. However their method suffers from the fact that 
adjustment of the pressures during iteration has to be done 
manually and required experienced guessing, and that it is 
only applicable for very heavily loaded contacts. One useful 
relation obtaining from this work in terms of dimensionless 
quantities 1s 
iy eSs0i298 5 (ate)? a S( US) TACwe ses? (1753) 

Osterle and Stephenson [8] obtained elastohydrodynamic 
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analysis included not only the effect of pressure-dependent 
viscosity and elastic distortion but also the effect of 
compressibility of the lubricant. However, no results were 
obtained for 'soft' elastic materials. 

Herrebrugh [9] showed that the hydrodynamic and 
elasticity equations in elastohydrodynamic lubrication can 
be combined into one single integral equation. A numerical 
solution of this integral equation for constant viscosity 
was presented. However, the resulting solutions were only 
claimed to be approximate. 

Cheng [10] presented the isothermal elastohydrodynamic 
theory for the full range of pressure-viscosity 
Pparameter,g*. His work suggested that for small values of 
g*, the lubricant could be treated as isoviscous fluid. 

In the late 1960'S, increaSing interest was given to 
the elastohydrodynamic problem lubrication for materials of 
low elastic modulus. Roberts and Tabor [11] employed an 
interferometric technique to study the fluid film 
lubrication of rubber. He showed that the deformed shape of 
rubber surface took the shape of a tilted pad surface along 
most of the contact region. 

Higginson [12] presented a model experiment in 
elastohydrodynamic lubrication of rubber. From his pressure 
measurements, he found the development of a near-Hertzian 
pressure curve as load increased. 

Roberts and Swales [13] studied the elastohydrodynamic 


lubrication of a rubber cylinder sliding over a glass plate 
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using an optical interference technique. They also used an 
iterative process similar to [7] to solve the basic 
lubrication and elastic equations simultaneously and 
obtained solutions mainly for lightly loaded contacts. They 
found the discrepancy between theory and experiment at light 
loads and high speeds, this might be accounted for by 
side-leakage effects. 

Baglin and Archard [14] presented an analytic solution 
of this problem, which was similar to the type of Grubin 
theory except they used simplied Hertzian displacement 
equation. In addition, they required to find a criterion for 
the pressure at the inlet edge of a nearly Hertzian region 
to ensure that the necessary distribution of pressure was 
generated within the Hertzian zone. Their results of the 
minimum film thickness had to be modified by a factor of 
0.8, in order to provide a comparison with the earlier 
computer solutions [9]. 

Swales, Dowson and Latham hy S:) described a 
comprehensive experimental and theoretical treatment of the 
hubsicationhof softtieelasticamaterials. cTheseffectptof meiniet 
Starvation on oil film thickness was also investigated. The 
experimental study involved a sliding contact between a 
Rubber cylinder .and a smooth plane. Pressure *and film 
thickness profiles were measured at the same time, using 
piezoelectic and capacitance transducers respectively. The 
theoretical work was similar to ‘that described on detail 


elsewhere. «[7]i, @#except ‘the lubricant was assumed to.be 
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isoviscous. Again the method suffers the same difficulties 
as stated previously. In addition, they found that it was 
difficult to achieve solutions for the case of large surface 
displacements. They were unable to decide whether the reason 
for the instability lay with the equations or with the 
computing methods. In general, they found both the 


experimental and theoretical results were comparable. 


1.3 Equivalent cylinders 

In the present study, a film of relatively short width 
will be considered between the contacts. Therefore the 
undeformed solids can be sufficiently represented by two 
cylinders in the region of the contact zone. These _ two 
cylinders can mere nes be represented by a geometrically 
equivalent cylinder near a rigid plane as shown in Fig(1). 
The geometrical requirement is that the gap between the 
cylinders in the original and equivalent contact should be 
the same at equal values of x. 

If the centres of the two cylinders lie on the same 
Side of the common tangent at the contact point then the 


radius of the equivalent cylinder takes the form of 
=—- - where R,>R;z 


Force components on the initial and equivalent systems shown 
in Fig(1) can be related as follows: 
For hydrodynamic force in vertical direction: 


Wy 1=Wy2=Wy=J p dx io) 
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hy,=R,-{Ri-x? = x?/2R, 
Similarly 

h2, = x?/2R2 

Hence 


h=h, +r) 1 +he 


=het(x?/2) (1/R,+1/R2) 


=h, +x*/ ZR 


where 1/R=1/R,+1/R2 


Fig(1) Representation of a contact between two 
cylinders by an equivalent cylinder near a plane 
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For hydrodynamic force in horizontal direction: 


Wie -) p-dhi=-Ci px ox eR: (1.6a) 
Wyx2=-J p dh2=-(J px dx)/Rz (1.6b) 
w,=-JAp. dh=<(fapx axe ZR C4664 
Hence 

Wx1=(R/R,)w, and wy2=(R/R2)w, (1.6d,e) 


For drag force components: 


F2=F =J tr dx=J ul(du/dy) | dx 61.7a) 
h h y=h 


Fy=Fio=s To ax=/ WUC W/ OF xa ax Ces) i) 


1.4 General assumptions 

The assumptions which will be made in the analysis are 
stated below 
(1) The cylinders are infinitely long i.e. side leakage is 
neglected and flow becomes two dimensional. 
(2) All inertia effects are neglected. 
(3) The width of the contact region is much smaller than the 
radii of the cylinders. 
(4) The film thickness is much smaller than the width of the 
GOnRTACT. 
(5) The flow is in steady state. 
(6) The lubricant is incompressible. 
(7) There is no slip at the boundaries. 
(8) The solids will be treated as perfectly elastic and ina 
condition of plane strain. 
(9) The contact will be represented by an equivalent 


cylinder near a plane. 
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2. PHYSICAL MODELLING APPROACH IN SOFT CONTACT 


2.1 Description of the elastohydrodynamic model 

A physical elastohydrodynamic lubrication model as 
shown in Fig(2) is postulated from which a numerical 
solution of finding the minimum film thickness is presented 
under isothermal condition. The theory of elastohydrodynamic 
lubrication of high elastic modulus | materials was not 
developed until Grubin [6] included the effect of elastic 
distortion of the bounding solids and the pressure-dependent 
viscosity of the lubricant. However, little attention has 
been paid to the development of the theory for _ soft 
materials in which deformation rather than 
pressure-viscosity effect becomes dominant. The lubricant in 
such case can be regarded as isoviscous fluid as suggested 
by Cheng [10]. 

Since there is no significant increase in pressure due 
to pressure-viscosity effect in the case of soft contacts, 
the generated pressure of the lubricant must be entirely due 
toy hydrodynamics action, Ay. ene, a converging bearing 
configuration, must -exist. along the pressure zone. In the 
light of thas-configuration ‘and the interferometric study 
done by Roberts and Tabor [11], it is postulated that the 
deformed cylinder under elastohydrodynamic lubrication will 
Have “a. Shapes like tilted-pad surface of finite ‘stepe along 


most of the contact region as shown in Fig(2a). 
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Fig(2) Geometry of lubrication film and 
model pressure distribution 
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Furthermore, from the film pressure measurements done 
by Higginson [12], pressure profiles for increasing load 
were found to be very close to that of Hertzian pressure 
distribution as shown in Fig(2b). In the ivght ‘of Pgthiis 
experiment, it is postulated that the deformed shape of the 
cylinder outside the tilted pad region will be the same as 


thet in a@dsy»contacti 


2.2 General equations for fluid 

The general equations governing the steady, laminar 
flow of an incompressible lubricant through an infinitely 
long rotating cylinder and a plane slider, as shown in 
Fig(2), are the momentum equations in the x and y directions 
and (aéthee contimuntyreequatrnenicifhes general -form of the 
equations may be listed as follows: 
x-momentum equation 

Oumu~oU- 1 op 62007 U 
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y-momentum equation with body force neglected 
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Continuity equation 
du OV 
—+—=0 C237) 
dx dy 
ipa) Sithe later Janalysis;,- the “following@velocity 


boundary conditions will be adopted. 
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y=h : u=-u2,v#-u;z(dh/dx) 6264c da) 
Particularly for the presently discussed model, the 


following pressure boundary conditions will be adopted. 


xetoo:. p=0 e205a)) 
dp 

x=0 : p=p =p, and —=0 G25C7,0) 
max dx 


2.3 Specialized equation for lubrication 
For representative lubricating films, the inertia 
forces can be shown to be negligible compared with the 
viscous force by the following order of magnitude analysis. 
Selecting appropriate reference quantities, the 
variables can be expressed as follows: 
x=xL, 
y=yh, 
u=uu, 
V=EVV, 
P=PP, 
Now substituting alleuthecstransforméd variables into the 
governing Equations(2.1) tove(2.3)4 and=.comparing the 
relative magnitude of the viscous terms with the inertia 
terms in each of the momentum equations, we obtain 
For continuity equation 
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be used in the order of magnitude analysis for the 


equations. 
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au/ax and avAdy areaof order’ 1: 


Wevmayen. must. “be of order 1. 


VpbezUeNe=l, we can Write v,= u,(h,./L,) which will 


momentum 


For the x-momentum equation 


du poum~oe op «1 bee Otro ou 
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Since the lubricating film is very thin as compared to 
the width of the contact region, i.e. h,/L,<<1, therefore 
the viscous term 0?u/dx? is negligibly small compared to 
0?7u/dy?. If Rex<<1, then the inertia terms can be neglected 


compared to viscous term. 


Hence only pressure and viscous 


forces will be present. Consequently Equation(2.1) reduces 
to 
Oprme 2u 
wy, Zino) 
Oxmady? 

In fact, Rex is much less than one in most’ lubrication 


Situations. 


For the y-momentum 
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Now comparing the terms on the right hand side of 
Bquations(2.6)@and (227);) we have 


From Equation(2.6) 


From Equation(2.7) 
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oy Ly He dy? 


Pe UiS Moted Sthat Goep/dy is -(h,/Ly) times the order of 
magni tudemootGa dp/erta Since, “hy/Ly<<1>e thecGvariation of 
pressure across the lubricant film is quite insignificant. 
Therefore the momentum equation for the y-direction can _ be 
omitted altogether. Thus write, dp/dx=dp/dx. 

Now the reduced set of equations only includes 
x-momentum equation and continuity equation of the following 
forms. 


x-momentum equation: 


dpapdGu 

Livitied (228) 
axGkoay? 

Continuity equation: 

du dv 

Sta: () (2.9) 
Ox Oy 


Furthermore, Equations(2.8) and (2.9) can be combined 
into one single equation as described below. 


Integrating Equation(2.8) twice with respect to y, yields 
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Applying boundary conditions(2.4a,c), leads to 
1 Pedp.vh 

c,=-(u,;-u2z)-—- — — and C2=-U; 
h yp taxG2 


Hence the velocity distribution is 


tiedp.y* hy y 

u=— —[—-—]+—(u,-u2)-u, (2510) 
PROX J2ma2 h 

Integrating Equation(2.9) with respect to y, leads to 

geeh oh h 

—J u dy-(-uz)—t+[v] =0 

ox 0 OX 0 

Applying boundary conditions(2.4b,d) the second and third 

terms of the above equation will cancel each other out. 


Substituting the expression for u in Equation(2.10) into the 


above equation, and then carrying the integration, leads to 
aa eee ne 


which is the famous Reynolds equation. 

Upon integrating Equation(2.11) with respect to x, leads to 
imaoenwsU 

mest eines —=-—ht+c 3 
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Now applying boundary condition(2.5d), leads to 
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where h, is defined as the film thickness at which pressure 


is maximum. Hence 
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2.4 Film thickness equation with Hertzian relations 
In region I, as shown in Fig(2a), the film thickness 


variation is linear with x and of slope m. It takes the form 


h=h,,; ,+m(xta) G25133 
Hence 
hy=hs;,*ma (2244) 


In region II, the film thickness will take the shape of 
Hertzian displacement. 
h=h,+6 
But from Equation(2.13) 
h;=hnint2ma 
Therefore 
h=hmint+2mat+d 22S) 


From Hertzian relations[16] 
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2.5 Non-dimensional procedure 
Using the Hertzian relations in Equations(2.16), the 
Equations(2. 12)% (2:13) and G2aeS9 can be 
non-dimensionalized by introducing the following 


dimensionless variables. 
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a R R E' E'R E'R 
Equation(2.12) is non-dimensionalized to 
dp* 8w* h*-h} 


=-6U*(——) 1/? 
ax; T eee 


) (2a) 


Substituting the dimensionless form of Equation(2.14) into 
the above equation, yields 
dp* 8w* We Shit aan C8 wt): 407 


=-6U*(——) 17? (———____________—__ ) (2 be) 
ax: T mt 3 


Equation(2.13) representing film thickness in region I is 
non-dimensionalized to 
8w* 
Pee a ptm (——). Sot Pees), ea a soir) C2a018;) 
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Equation(2.15) representing film thickness in region II is 


non-dimensionalized to 


Sw* 
h*=h,,; ,+2m(——) '72+8* , Cx sat) CZrsic)) 
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The corresponding non-dimensional boundary conditions become 


Xtato p*=0 (2.20a) 
= Ah «hep TS0 (2.20b) 
x ;=0 3 p*=p* =P6 G2A20c) 
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2.6 Numerical analysis 

In this section numerical methods will be described for 
the solutions of hades and m through the pressure 
equation(2.17b), subject to the boundary conditions(2.20), 
where the film thickness is given by Equations(2.18) and 
(2.19) in regions I and II, respectively. Input parameters 
U* and w* must be given to define a particular problem. 

Initially, the pressure equation is solved in the first 
half of region I by first assuming hj,;, and m, and applying 
boundary conditions(2.20b) at xi=-1. The pressure equation 
.1is then integrated, using Simpson's rule, in the positive 
Me-avrect tonmand boundary -condition(2.20c) isicheckéd. Ff it 
is not satisfied, a new value of hjfj;, is tried, but with m 
remaining constant, until the boundary -condition(2.20c) is 
Satisfied. Then with these values of new hf,, and m, the 
pressure equation(2.17b) is again integrated in the positive 
x=direction between xees J and X,=to, and boundary 
condition(2.20a)*%isichecked? If it is not satisfied, a’ new 
value of m is tried. The whole process is repeated until the 
assumedivalueswor h,;, and m Satasfy ‘all. tthe boundary 
conditions (22.20). 

In actual numerical procedure, it is assumed that p”=0 
at xt=t+t» reaches its asymptotic value when the magnitudes of 
the pressure and pressure gradient become less than a 
Specified sfraction .(10°°)) of: the maximum pressure. and 
maximum pressure gradient, respectively. Also the relative 


error allowance on matching p*(max)=po at xj=0 is specified 
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to be less than 107°. 

A Newton Raphson iterative technique is employed to 
converge hyin and m on the solutions. This technique, which 
1s quadratically convergent, can be used because of the fact 
that a pair of nonlinear equations can be set up as the 
functions of two variables, hfi,, and m. Basically, this pair 
of nonlinear equations iS derived from the boundary 


@eanditions(2.20) and 1s of the following form 


x7=0 dp* 
|B al 9 ee Rees 110 | [ Weds n=, (22a) 
xnSs-dredxe 
uf 
and 
xi=1 : Xo=to dp: 
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PEGES noted that in order to satisfy all.i the "boundary 
conditions(2.20), the assumed values of hj;, and m must be 
gine roots of the “functyrons f, ‘and: f2, 
teen le eens ast) =bo (hae me Os 
The roots are found by Newton Raphson iterative method as 


described below: 
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All the partial derivatives are approximated by making 
a small change in the value of a variable and dividing the 
change in value of the function by the change in value of 


the variable. Mathematically, it takes the form 


ah} 13 BA fe. 2 tm aphege teige: tm) hey i 
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The subscript t indicates the number of the iteration. 
The iterations are carried out until the difference between 
two consecutive sets of values is less than a_e specified 
relative error allowance (107°). 

A criterion for ‘convergence of Newton's method 
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on an interval about the root, the method will converge for 
any initial values of hj, and mein the interval. 

The method of choosing the initial guesses of hj;, and 
m is based on the physical insights from Equation(2.17a). We 
observed that decreasing chfin,owith’mcheld constant, will 
increase the entire pressure distribution. It is because the 
term (i7hat) eink Equation(2.17)» becomes large as hnin 


decreases. On the other hand, decreasing m, with h,;, held 


hs 


= 


prise yd betamizozage- ete eau hte iad Seba gen 
= 
nay petbevib bes sideiver 6 te selisv- she mins 


+ Bs 
i906 ©6spwisvy )6fnt)6 6 .spmaeno end ve noi sony eit 310 4 can 
ae. 


vs 


miot ef4 geNsey te ,ylipeisensdaelt 
It (wie 


a : (heh 4 eel li wit ame 
Se a a 


ue 


re ; 
ia. rach ; visoibda $s. syeisedes set 
a. aes 223 5h eat iy. 2uS bs. xTSr 53'S ens i 8 2eRE 
229) ei eau lev Yoretseaz ow! Sasi 


a4 SP20.80K! 5 Tate . 
| ay 


2'naves\ +9 S>nepteUloo. (aes te Pesare ae 


7) +nAs BS) Yin noid rSnaa a of 


rf 
1 


a ae 


ied agsevaeg Lite Sodzem aes pei © sii4 Sad lavas I Re wu 
Leviedin Sha att. _ roam 20! Sand aw ee y 


a 
7 ” 
ane ,.eat to 262 3up: fehain i Siz davandis 3 bos isa eae 7 


aW .(e80l. S$) cd teupe me3t) aero! ated, fapteyda sid no sean: & 


_ 


* ci 7 
(Liw toatenes Slade a8iw , , aH oiiasersae tsdz te — 


SAY aevered @r 1 .nolsudii yeets S3udestG) avlins ods 
sind 28 spvel  gsmoosd (Pr Aline t Seip ab. Le dhe) 
~~ - 7 


Slee inl datiw on babzasyss5 , brad. yehia sn3. a0. Jame 


ez 


constant, will lower the entire pressure distribution. It is 
because the term (h*-hz) in Equation(2.17a) becomes eaait as 
m decreases. It was found that pressure was more sensitive 
to Nhmin than to m. After few trials, better initial guesses 
Satisfying the convergence criteria(2.22) are obtained and 


then used in Newton's method for finding convergent 


solutions® ofihti eeands me 


2.7 Results and discussion 

Results obtained from the previously described method 
are presented in Figs(3) and (7). Input data given in 
Appendix I have been used to obtain these results. 

ReferrangioterFigs(siethroughsd5) ,aitedisenoted that) at 
a given speed, the minimum film thicknesses for both rigid 
and elastic cylinders are decreased as load increases. At 
very light loads, predictions of minimum film thickness 
obtained by rigid cylinder theory are in good agreement with 
the present results. It is because surface deformation under 
such condition 1S insignificantly small as compared to the 
thickness of the film. However, as load increases, 
deformation becomes increasingly important. Under such 
cases, the minimum film. thickness obtained from the 
elastohydrodynamic model is found to be much higher than 
thatepredretecm aby. rvgqid @cylinder -theony. “lies "simply 
because cylinder is deformed at increasing loads. It is also 
observed that ,at higher loads, the minimum film _ thickness 


becomes very insensitive to change in load. 
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It 1S indicated, from Fig(6), that at a given load, the 
minimum film thickness is increased as speed increases. This 
increase in minimum film thickness due to change in speed is 
found to be quite uniform along any given load. 

In order to examine the accuracy of the present 
results, it is necessary to compare the presently obtained 
results with the experimental data [15]. For the sake of 
Conpannconenncye7 (U2 “re plotted against vwi/(U')ia* “as 
shown in Fig(7). The discrepancy between the present and 
experimental results at light loads and high speeds is due 
to the effect of side-leakage. Whereas, the discrepancy at 
high loads reveals that the assumed film profile is no 
longer true. However, reasonably good agreement is found in 
the range of 0.9<w*/(U*)'’?<20.0. Most of the engineering 
applications in practice fall within this range of loading. 
Therefore, the described model, based on the concept of 
Hertzian deformation associated with a tilted-pad 
configuration shown in Fig(2), is proved to be useful in the 
analysis of the elastohydrodynamic lubrication. of. soft 


materials. 
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3. THEORETICAL SOLUTION IN LAMINAR REGIME 


3.1 General approach 

A theoretical solution to the elastohydrodynamic 
problem must consider the basic equations of elasticity for 
the solid, and equations of momentum and continuity for 
Eluid. Though the analysis in [Chapter 2 yields a useful 
approximate solution for film thickness and pressure 
distribution; it does not tell us the actual pressure 
profile and film thickness in the contact region. Thus, in 
order to examine the actual details of the problem in a 
variety of conditions, the basic elasticity, momentum and 
continuity equations must be solved simultaneously. The 
initial geometry of the physical problem as shown in Fig(8) 


will be considered. 


3.2 Equations for fluid 
As described in Chapter 2, the equations for fluid 
comprising of momentum and continuity equations may be 


combined to form one single Equation(2.11) of the form 


where U=u,+tu2 Caen b) 
This equation is based on the assumptions that the flow is 
laminar, l.e., sat low. Reynolds “number, the viSeOsuty 1s 
constant across the film, and the width of the contact 


region is small compared to the length of the cylinder so 
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that side leakage can be neglected. 
The above equation is subject to the following boundary 


conditions: 


dp 

X=Xo $ p=—=0 GS2ayD) 
ax 

Xero : p=0 Vae2o 


Equation(3.1a) is integrated with respect to x, to yield 


fadp hh U 
aj. ee Se ESS —ht+c, 


pHaxa 12 2 
Now applying boundary condition(3.2b), leads to 
U 


Cy=—No 
Z 


Hence integrated form of Equation(3.1a) becomes 
dp hoo u 


—=6Upn ( 
ax he 


) (earees)) 


At heavily loaded contact, especially in metallic 
eontact, pressure-viscosity effects Wadd arise. The 
viscosity may be considered to increase exponentially above 
their atmospheric pressure values, and of the form 
u=uoexp(ap) | (354) 
Introducing Equation(3.4) into Equation(3.3), yields 

| dp ho-h 


exp(-ap)—=6UUuo ( 
ax igi 


) 37055) 


Since one of the boundary conditions(3.2) has been used in 
fanding tthevconstanty cq, ABquationd3 5)\ae1sythenefore . only 
subject to the remaining two boundary conditions: 


MaKe: p=0 (3.6a) 
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x=to : p=0 (3.6b) 
Now it 1S convenient to introduce a new variable q [7], 
known as reduced pressure, such that 

dq he-h 


—=6Uuo ( 
dx ha 


) Gols) 


where 


dg dp 

—=exp(-ap)— (3.8) 
dx ax 

Integrating Equation(3.8) with respect to x from the outlet 
point ,Xo, to some general point ,x, we obtain 


g Pp 
f “dg=) exp(-ep)dp G3229)) 
Gxo Pxo 


Fromenthes boundany® condition(3u6alieipgs tsizero.aréi qj¢ is 
set equal to zero, the function q will satisfy the same 
boundary conditions at xX, as the pressure, p. Thus from 
Equation(3.9), we get 

1 
q=—[1-exp(-ap) ] (3.310) 

a 

It can be readily seen that Equation(3.10) also 

Satisfies the second boundary condition(3.6b). Hence the 
function q satisfies the same boundary condition as_ the 
pressure p, that is 
ya cee Ue Ho 10) (a3 aba) 
x=to ; g=0 (373169 
Thus this problem can be solved independently. Hence, once 


Equation(3.7), respresenting the pressure distribution that 


would*OccUrein anlisoviscousyilubricanteeofiheviscosity o,1is 
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solved then the actual pressure distribution can be found by 
Equation(3.10). That is 
1 1 


=—ln( 
a L=ag 


) O32) 


Eoveciont3612)- (is Svalid only: when aq<1, i.e. g<1/e: If 
q>1/a, the actual pressure goes to infinity. Physically, 
pressure will never become infinity because high pressure 
will deform the surfaces into Such a shape that the pressure 


remains finite. 


3.3 Elasticity equation 

in, iebastohydrodynamic “lubrication “situation, it is 
reasonable to assume that the bounding solids are " a 
CONdition Of plane’ Straims It 1s due to the fact that the 
width of the pressure region is very small compared with the 
radius and length of the bounding solids. So the stresses 
and displacements are uniform along the length, except near 
the ends. In addition, this fact further allows us to assume 
that the surface displacements and stresses of the bounding 
solids can be approximated by those of semi-infinite solids 
subject to the same normal load. The final film thickness is 
therefore obtained by adding the vertical displacement to 
the initial thickness of the film. 

First =’consider"a line load on“a"*semi-inmftinite “solid as 
shown in Fig(9). 

Under «plane ystraim.’) conditzon, 1,¢€. w=constance and 


a¢/az=0, where $ is a dummy variable, we obtain 
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Bor drrect strain 
ou OV OW 
Ce Saray Gn oe ee ’ €,=—=0 
OX oy OZ 
For shear strain 
OU OV du OW OV OW 
eye te ya he Oe ieee es —=0 Pot veen —=0 
dy Ox OZ 20% dz Oy 
Therefore only stresses o,,0,,0, and_r,, are required to be 


considered. 


Hooke's law relates strains to stresses as follows: 


1 dU 

€,=-[o,-u(o,to,) J=— (25 13a) 
E Ox 
1 OV 

ey=-[Lo,-u(o,t+o0,) J=— eRe) 99 
E Oy 
1 Ow 

e,=-[o0,-v(o,+0,) ]=—=0 C3.n1Scy) 
E OZ 
GR ES (wer ve) OU OV 

Fe te ees (3.13) 
24g E dy ox 


From Equation(3.13c) yields 

Te viG +O) | (3,14) 
For a solution tottexiispaAmtebothtsequilrbrium and 

compatibility equations subject to the appropriate boundary 

conditions must be satisfied. 

Equilibrium equations in two dimensional form are 
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COUT s F 
ar =Q (3.15b) 
ey s0d0x 


Compatibility equation in two dimensional form is 


Orted? 
+ 


axgnay? 


( )(a,+0,)=0 3eoiG)) 


It is convenient to introduce a stress function ¢ (x,y) 


such that 
07 076 07@ 

C= pe Oy pe gg ee Cee er: 
aye Ox? Oxdy 

Equations(3.17) automatically Satisfies equilibrium 


equations(3.15), but the eompatibility.! Equation(3.16) 
becomes 
0*¢ 0*¢ 0*¢ 


+2——__—__+——-= 0 03 998) 
On. SOKROY 8a y* 


The following Boussinesq stress function satisfies the 
equilibrium and compatibility equations for the case _ shown 
nei g:( 9). 


wr @sin @ 


Then we have the following stress components 
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These stress equations satisfy the natural boundary 
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conditions, namely stress free at the surface, i.e. one at 
6=90, except at r=0, and 0,70 as rIm, 

It is convenient in the present problem to use 
carteSian coordinates for the stress function and stresses. 
Stress function becomes 

W x 
=-—x tan’! — 
us yi 


and stress field 


07d 2wx’y 
05 = ca eae ate are (3194) 
oy? mike? yd) Z 
OfOS L2WYE 
Oy= LS OCT ET er (3. 19b) 
0x2 TARRY Ce 
-07¢ 2wxy? 
Tx y= Bese (3%19¢)) 
ax ey kim (7 bye) 2 
Introducing Equations(3.14) and G3 19) into 
Equation(3.13a,b) and integrating we get 
Ww x x 
u=-—[(A+B)tan™'! — - (A-B) ]+£ Cy) (3% 20a) 
T y Mew y7 
W : x? 
V==—fAiIn(xet+y2 pass} 1#6(x) (2°.20b) 
1 KcHy < x 2hy:* 


wheres Etylabandieilalemanés fencieons of y )and ix ‘only, 
Rogue: sacae ands A=(1-v2)/E, B=-utii1+u)/5. | 
Now by symmetry at any value of y 

-u(-x)=u(x) 


Hence £(y)=0 
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Tocfind-£(x), iBquations(3 713d) and (3.19¢} 
du OV E 2wxy? 
Te) ore ee oe bh) 
Differentiating Equations(3.20) and then introducing them 
into Equation(3.21), yields 
d 
—f(x)=0 
ax 
Hence f(x)=constant 
Therefore, Equations(3.20) give the displacements anywhere 
an eenesbodyrdue to a cline Loadv:at the origin. But: we are 
primarily concerned with the vertical surface displacement, 
WrevVvat y-=0,--due to alsurtace eee tre p(s) between x=s, 
and x=sz as shown in Fig(10). This displacement can be found 
Byereplacing- wala xivin (3.205) 1 with (p ds6))) and . (x-s) 
respectively, and then integrating with respect to s, we get 
(107) ais: 
v | =-————. {___ p(s)1n(x-s)? ds + constant (3.22) 
y=0 TE S 
In the contact between two solids, both surfaces will 
be deformed by the surface pressure, and if the subscripts 1 


and 2 are used to denote the two contacting solids, the sum 


of the surface displacements become 


4 S2 

v(x)=-—— J Di sj im(x=s) ds) saconstant S23) 
Te S14 

where 

1 1 St 48) (te 52) 

SS | tLe 

| 2 Ey E 2 


Elasticity equation similar to the one derived above was 
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also used by Dowson and Whitaker [7] in their analysis of 
surface displacement. 

Pom wee noted. ‘that, sfrom Chapter “1; Equation(1.4) 
represents the film thickness in lubricating rigid contact. 
Therefore, the resulting film thickness in lubricating 
elastic contact can be obtained by adding the surface 
displacement to the undeformed parabolic film thickness, 
that is 

x 
h(x)=h + — + v(x) (3.24) 
Eg 2R 
where v(x) is defined by Equation(3.23) and the central film 
thickness, he, can be determined from the lubrication of 


bigid contact. 


3.4 Non-dimensional procedure 
Babatilons (3 oe rian tie ets ke) and » (8824). “can. be 
non-dimensionalized by introducing the following 
dimensionless variables. 
xt=x/R, gt =q/Bi; po=p/E')ohe=ho/R 
hg=he 7k; h?=hYR, So) Sabu xGexeyR 
Vv =V/R, alee7R U*=Un"o/E'R 
Equation(3.7) is non-dimensionalized to 


dq* h3-h* 
i é ) (32.25) 
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which is subject to the boundary conditions(3.11) in 


dimensionless form 


+ 


X 


=to 


=x : q*=0 3026) 
q*=0 G3i7n2iGb) 


Equation(3.12) is non-dimensionalized to 


1 1 
i) (oz aN) 
loc 


ee : 


g 


Equation(3.24) is non-dimensionalized to 


ar 
h* (x*)=hi+——+v* (x") (3. 28a) 
; 
where 
4 
Vax )=——2p xs yt 6c (3.28b) 
T 
and 
S2 
BGxsy yas ip (ste) ini seers oi ass ote OG) 
Sj 


The constant, C, may be determined by specifying the 
minimum film thickness as will be described in the next 
section of numerical method. Therefore, from Equations(3.25) 
to (3.28), the solution depends on three parameters, namely 


oie; BUR tand eh iw 


3.5 Numerical Analysis 

Numerical techniques of finding the solution. for the 
pressure equation(3.25) Subject to the boundary 
conditions(3.26) with the film thickness equation given by 


equations(3.28) will be described. 
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3.5.1 Solution for the pressure equation 

It 1s noted that from the boundary condition(3.26a), 
Phessoutleteupornt; }.x5; .isGunknown. The Waynto locate sxs uis 
that by first assuming a value of x$ and applying boundary 
condition(3.26a) at that point when the pressure 
equation(3.25) is integrated in the positive x*-direction. 
Botnda dy condi téon (3e26b)m isicthensz-cheekedy If it is not 
Satisfied, a new guess of xj is tried. 

The method of choosing and correcting the assumed value 
of x6 is based on the fact that moving xg to the right will 
decrease the entire pressure distribution, whereas moving x3 
to the left will increase the pressure distribution. This 
will allow us to find upper and lower bounds for xj. Since 
if the assumed value of x$ is too far right to the actual 
value, then pressure becomes negative while integrating in 
the poSitive x-direction; hence an upper bound. On _ the 
contrary, if the assumed value of xj is too far left to the 
actual solution, then pressure becomes too high at_ the 
inlet; hence an lower bound. Once the bounds for outlet 
point are found, the solution for xg can be easily obtained 
by a simple bisection method. 

In actual numerical procedure, it is adopted that q*=0 
at x*’=+t» reaches its asymptotic value when the magnitudes of 
the pressure and pressure gradient become less than a 
specified fraction(10°*) of the maximum pressure and maximum 
pressure gradient, respectively. Equation(1.1a) in Chapter 1 


will give a rough guess for the outlet Dotmertx>o: 
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Te should bee enoted| that’ the “f2lm-.thickness is 
considered to be fixed when xj is being iterated with the 
pressure equatiron(3.25)2 Once x5 is found, the actual 


pressure distribution, D3 can be obtained by 


Equation(3.27). 


3.5.2 Solution for the film thickness equation 

Once the pressure distribution is determined for a 
given initial film thickness, it is timé'to evaluate a 
corresponding film thickness variation from Equation(3.28a). 
This requires the evaluation of the deformation integral 
Dex Je yin ‘Bahavron¢se2oe. ssi hevexpression,;“o(x ) + contains a 
singularity -anithe antegrand, f.e7 © In @’-sit)o-e vas» -x*35", 
and so it cannot be evaluated by straight forward numerical 
integration. This difficulty is removed if the pressure is 
expressed by a function which enables the integral to be 
performed analytically. In the present study, D(x*) is 
evaluated by dividing the pressure curve into Suitable 
segments. Each segment contains three pressure points and is 
represented by a second degree polynomial of ene form 
pris )eatbs +65. (3.29) 
where a, b, and cC are constants for one segment. For the 
purpose of computation, the origin 1s taken "at the second 
Pressure “point of each segment, as shown in big (id) @. Three 
pressure points of known positions are allocated “to each 
segment and their respective coordinates referring to the 


new origin can be easily found. With respect to the new 
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Origin, the new coordinates for the three pressure points 
beeoment (so, ppeakeurte.p; hivamdt i6s3 ap} sh). eaTheypoiare » ‘then 
substituted into Equation(3.29). Consequently, we obtain 
three similtaneous equations from which constants a, bandc 


can be determined in each segment. The constants are found 


to be 
@=pj 
Piianatp 2t7a | 
C= : = : 3 : 
S2 S; Sa Sy 
Dna a 
= = Sra KS 
S2 
With the pressures, piulsth in the Porn of 


Equation(3.29), we can now integrate Equation(3.28c) without 
any difficulty. The surface displacements due to each 
individual pressure segment can be found as described below: 


For each segment, we have 


So : 
Dxe as =| atbsetose 7) in (cies) ds" 
Sj 
whereonx’ =X 3-274) -s#=xi'. {Sx Pland siex | arex) 
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pressure segment. Therefore, the actual displacement due to 
the entire pressure distribution at the same point x* can 
then be determined by superposing the results of D(x*) 


obtained from each pressure segment, i.e. 
DCs?) (S230) 


where kK = total number of pressure segement. Rearranging 
Equations(3.28), we get 


xnee4 
== px" )4 
T 


Gr Cen + 
ee 


Using Equation(3.30), the function (h*-C) may be evaluated 
at every partition point. Then, after the minimum value of 
this function is found, the proper constant C may be added 
to the function so as to give the desired minimum film 


thickness, hjjn. 


3.5.3 Film thickness iteration 

With the numerical methods described in the above two 
SECCIONS 4103p al) pandu. (shovel ewe tare able “torfind=the 
pressure distribution from a given film” thickness and the 
film thickness from a given pressure distribution. This 
procedure allows us to check any given film thickness to see 
if it is a compatible elastohydrodynamic solution. That is, 
if a film thickness hj is given, the corresponding pressure 
distribution can be obtained from Equations(3.25) and 
(3227). The pressure is then used to “generate a new= film 


thickness h3 through Equations(3.28). If hj equal to h3, 
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then h; and its associated pressure distribution satisfy the 
elastohydrodynamic equations. If h} does not equal to h3, a 
method of choosing and correcting h} must be developed so 
that compatible solution may be obtained. In the present 
analysis, the resulting film h3 is used as the new guess for 
hj. This scheme of iterative procedure for film thickness is 
basedvnon tthetstagctetthate:ifiethe pressure’ is tootvhigqhtatra 
given point, there will be too large a deformation and _ the 
film thickness will be larger. Thicker films usually result 
in lower pressures. Hence there is a self-correcting feature 
for this iterative scheme which should lead to a convergent 
solution. It was found that the above scheme did not bring 
convergence in the cases where deformations were large 
compared to the film thickness. The problem of divergence is 
removed by introducing a weighting factor, 6, which weighs 
each new film with the previous film. This will prevent too 
large a change in film thickness in one iteration. Now 


(hj) =ph>t (1-8) Che) Co 3is) 
new one 


Where thee welgntinoqe tactor ps1," end  ~=4 srepresencs: che 
original iterative scheme. 

As mentioned ‘before in “section (3.5.2). the proper 
constant C may also be added to the new hj in Equation(3.31) 
to make the minimum value of hj agree with the desired 
minimum fam thickness, hjin, before a new pressure 
distributvon 1S _being iterated. 

It was found that either the film thickness, or the 
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distribution corresponding to a given film thickness can be 
used to check for convergence. 

The rigid «film profile can first be employed to 
initiate the iterative procedure. However, several 
iterations are required before a convergent solution can be 
obtained. 

Once the convergent solutions of the film thickness and 
the corresponding pressure distribution are determined, the 
load capacity, velocity field and drag force can be obtained 
from Equations) hh1t5) ¢inbeed0)weand:(1.7),, respectively. The 
dimensionless form of these equations can be written as 
follows: 


Dimensionless load capacity is of the form: 
wi=f p*dx* (Grs2)) 


Dimensionless velocity is of the form 


1 dp* ve 
u*=— Cyt eho et 
Pash ob ail he 


Gua) —U4 (eldec e's 


Dimensionless drag is of the form 


+00 Ow 


ES) axe (3.34) 


3.6 Results and discussion 
Ihesthis «section; eresults, obtained. from: theapreviously 
described method will be presented for both low and _ high 


elastic modulus materials under four different conditions, 
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(in) rigid cylinder with constant viscosity. 

(11) rigid cylinder with pressure-dependent viscosity. 
(iii) elastic cylinder with constant viscosity. 

(iv) elastic cylinder with pressure-dependent viscosity. 
Input data given in Appendix II has been employed to obtain 
these results. 

As mentioned in the previous section (3.4), results of 
film thickness and pressure distribution are dependent upon 
three parameters, namely 
(aii) minimum film thickness parameter, hjij,. 

(ii) speed parameter, U*. 
(iii) pressure-viscosity parameter, g* 
In the present study, the influence of hy,;, and g* upon 


film shape and pressure distribution will be examined. 


3.6.1 Results and discussion for soft material 

For materials of low elastic modulus, i.e. g* becomes 
small, Figs(12) through (14) show some pressure 
distributions under four different conditions for three 
different minimum film thicknesses. To examine the 
Significance of the pressure-viscosity effect, the pressure 
distributions obtained from constant viscosity are compared 
with that from pressure-dependent viscosity. It is noted 
that there rssnoedifference win pressure’ in both ‘cases, This 
means that the pressure-viscosity effect may be neglected 
and only the deformation effect need to be considered in the 


analysis of hibmrcatn om fo@tscoity contact ominence;° the 
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lubricant in* such case can be treated as isoviscous fluid. 
This supports the argument made previously in Chapter 2 that 
fluid is isoviscous. It is also noted that the pressure 
distribution obtained from elastohydrodynamic theory is 
found to be higher than that predicted by rigid hydrodynamic 
theory. The differences are mostly prominent at small values 
SCH Vie. 

It is observed from Fig(15) that as h,j;, is lowered, 
the, outlets pormtwstarts shiivingetocthe fleftaAbresult itiguin 
raising the entire pressure distribution. This is caused by 
deformation of the bounding solid due to the hydrodynamic 
pressure generated in the fluid film. Maximum pressure is 
also observed to shift toward the outlet as h,,, decreases. 
In addition, Fig(15) also shows that small decrease in low 
values of hpin can cause a large increase in the entire 
pressure distribution. 

Figs(16) through (18) show the corresponding 
elastohydrodynamic film thickness for three different values 
of hii» Fig(16) represents the one at lightly loaded 
conditions, where deformation is relatively small compared 
to the film thickness. Fig(17) represents the one at 
intermediate jodae) where deformation starts becoming 
important Macht isis -calsoghnotedesfromaurig(17)athat abndanrby 
tilted-pad region is apparent. This supports the argument 
made previously in Chapter 2 that tilted-pad region does 
exist for largersloads. Fig(18) represents’ the one atypsheavy 


loads, where deformation is relatively large compared to the 
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film thickness. The location of hii, is found to be shifting 
to the left as hai, decreases. This accounts for the shift 
of maximum pressure as mentioned in the previous paragraph. 

Velocity profiles at three different locations are 
shown in Figs(19) through (21). Referring to Fig(19), a 
linear velocity distribution is noted at the outlet. It is 
because the pressure gradient at the outlet is zero. In 
Fig(20), a bulgy parabolic velocity profile is observed at 
hn k« This shows the presence of favorable pressure 
gradient. Fig(21) represents the velocity distribution at 
whemeinlets @Back-tilow asenotedbatey?>ic2x107 *e#iIltcisiduesto 
the presence of adverse pressure gradient. 

Figs(22) and (23) respectively show the variations of 
load and drag for different values of hnji,. AS hna;jn lowers 
from 6x10°* to 4.5x10°*, the percentage increases in load 
and dragiarentoundytesbe. 752%eeapproximately ul But ~asaah, py 
lowers from 4.5x10°* to 3x10°*, the percentage increases are 
about 150%¢2sThusy avtrend ofsrapid rate offtancrease inveload 
and drag is noted as hjin gets smaller. The drag at y*=0 is 


found to be three times the drag at y*=h’. 


3.6.2 Results and discussion for hard material 


+ 


For materials of high elastic modulus, i.e. g becomes 
large, pressure distributions under four different 
conditions for three different minimum film thicknesses are 


presented in Figs(24), (25) and (26). To study the impact of 


thes Spressuresviscosity veffectssethe comparison between 
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pressure distributions obtained from pressure-dependent 
viscosity ‘and constant viscosity is made in the case of 
rigid cylinder. It is found that there is a significant 
increase in pressure due to pressSure-viscosity effect 
between x*=0.0 and x*=0.02 for all three values of hj,,. 
This effect is observed to be increasingly significant with 
Gecueasingi WininoeeOneethe ‘other— -hand; the effect of 
deformation with constant viscosity 1S examined. It is noted 
that pressure distributions obtained from deformation alone 
are found to be much less than that from pressure-viscosity 
alone. This shows that the pressure-viscosity effect is more 
influential than the deformation effect. 

Fig(24) suggests that for relatively high values of 
Data? thes effiece of fdeformapionycan be neglected. «Only sthe 
effect of pressure-viscosity iS important. In such cases, 
the elastohydrodynamic pressure distribution can be 
accurately obtained by considering the case of rigid 
cylinder with viscosity varies with pressure. However, 
Figs(25) and (26) show that as hjin gets lower, the pressure 
distribution, obtained by considering pressure-depentent 
viscosity alone, starts deviating from the true solution. 
Thisseshows ithat; seven) (though 9sdeformation is small, the 
effects of pressure-viscosity and elastic deformation are 
now becoming equally important. Therefore, in order to get 
ane ebastohydrodymamic jesohution, @uboth effects must be 
considered. The actual maximum pressure is found to be two 


times higher than that in the case of rigid cylinder with 
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pressure-dependent viscosity, and is four times higher’ than 
Ghatevin brigqidteandecelastic (cylindersnebothswithiconstant 
viscosity. This Significantly increase in pressure 
undoubtedly demonstrates the importance of 
elastohydrodynamic theory. 

Referring to Fig(27), a small decrease in low values of 
ay ecan give rise to a large increase in the entire 
pressure distribution. This increase in pressure, however, 
is not only due to the effect of deformation as described in 
the previous section EOG soft materials, but also 
Significantly due to the pressure-viscosity effect. Such an 
increase in pressure is found to be much drastic in the case 


+ 


OEthardematéerialymi..é€,=%q becomes large. The reason is 
Simply due to the influence of pressure-dependent viscosity 
in hard contact. It 1s observed that the maximum pressure 
and the outlet point are slightly shifting to the left with 
decreasing hy) 7. Butethiseshift ts notcasssharpias | othatviin 
soft material. The reason is due to the fact that hard 
material has high elastic modulus, i.e. higher resistance to 
deformation; therefore, the area subject to deformation is 
relatively small as compared to that in soft contact. This 
canibe#well@vilustrated an. Eigs(28) tomi3s0), 

Figs(31) through (33) respectively present the velocity 
profiles at three different locations. Referring to Fig(31), 
thesvelocity dietributtonhatithe outletr is#punimeugmatacrogs 


the fluid film, since there exists zero pressure gradient. 


Fig(32) shows that, due to the presence of favorable 
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pressure gradient, the bullet-shape velocity profile is 
observed at the place where film thickness is minimum. The 
maximum velocity occurs at the middle of the minimum 
thickness of the film and is about 1.6 times higher than the 
Sliding speeds. Furthermore, in order to satisfy the flow 
continuity, this local maximum velocity must be the global 
maximum velocity in the entire fluid field. Fig(33) shows 
the velocity profile at the inlet of the flow. Reverse flow 
1s observed at the central part of the film because of the 
presence of adverse pressure gradient. 

Referring to Figs(34) and (35), variations of load and 
drag .for different minimum film thicknesses are noted. As 
Kietiw LOWE ELSo Tt rOmass sox On eat On OL; the percentage 
increases in Loagdeaancemarag are about 39% and - 17%, 
respectively. Whereas, ae hear lowers. trom ax 107 > tomu6e5r24 
the percentage increases are about 128% and’ 735%, 
respectively. A tremendous rate of increase in load with 
smalk «decrease of low values of Aig is observed. The rate 
of increase in load is much faster than that found in soft 
material. The rate of increase in drag is Similar to that 
found (in sottwecontact. Liisi noted that thesdrag atmy =0. and 
y*=h* are equal because of the flow symmetry. 

Inforder to, study <the: ‘accuracy .of “the theoretical 
results, a comparison is made between the present 
theoretical results with the earlier experimental and 
theoretical results obtained by others. For the sake of 
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shown in Figs(36) and (37). Reasonably good agreement is 


found. 


4, CONCLUSION 

A study was made Le an isothermal elastohydrodynamic 
lubrication off materialsssof both» low’ and ‘high elastic 
moduli. The described physical model, based on the concept 
of Hertzian deformation associated with a tilted-pad 
surface, for materials of low elastic modulus was proved to 
be useful in predicting the minimum film thickness, 
Re jee movers “auwace (range of loads, wi/(U° !'’2 = 90.9 
to 20. In the range of conditions where elastic deformation 
was important, it was found that for a given load and speed, 
the minimum film thickness was much greater than the rigid 
cylinder theory would indicate. The increase in minimum film 
thickness due to change in speed was found to be quite 
uniform along any given load. 

Results obtained from the theoretical analysis 
indicated ~*that-7 fom Gubrrcation of soft materials with 
g*=0.933, the effect of pressure-dependent viscosity could 
be neglected. Only the effect of elastic deformation should 
be considered in such cases. However, for lubrication of 
hard materials with g*=3000, both effects of 
pressure-viscosity and elastic deformation were found to be 
equally important in the generation of the resulting 
pressure distribution. Over the range where deformation 
became important, it was noted that for a given minimum film 
thickness and speed, load capacity was much higher than that 
obtained by rigid cylinder theory. Drag and load were found 


to be increased with decreasing minimum film thickness. 
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Small decrease in low values of the minimum film thickness 
would cause a large increase in load. In other words, when 
deformation is significant, minimum film thickness was 


slightly dependent on load. 
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Fig( 34) Variation of load with minimum film thickness 
for elastic cylinder 
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Fig(35) Variation of drag with minimum film thickness 
for elastic cylinder 
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APPENDIX I 


Input data 


Radius of cylinder (R2) 

ple eT of cylinder 

Atmospheric een tare viscosity (no) 
Modulus of elasticity (Ez) 
Poisson's ratio (v2) 

Cylinder rotating speed (uz) 

Slider speed (u,) 


Load (w) 
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0.0635 m 

rubber 

0.138 Ns/m? 

24x10° N/m? 

oS 

0 

0. 598mse- 2.37 m/s 


525EN7m - 57800 N/m 


(em 


(ei need ani tsto® 


APPENDIX II 


Input data for soft material 


Radius of cylinder (R2) 0.0635 m 
Material of cylinder rubber 
Atmospheric absolute viscosity (uo) 0.138 Ns/m2 
Modulus of elasticity (E>) 24x10° N/m? 
Povsson'-s ratvo (v3) 0.5 

Cylinder rotating speed (uz) 0 

Slider speed (u,) P85 imis 
Pressure viscosity exponent (a) 455105 men 
Minimum film thickness (h,; ») 1.905x 1072 mm to 


BOL LO” “mm 
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Input data for hard material 


Radius of cylinders (R,, R2) 0.0508 m 
Material of cylinders steel 
Atmospheric absolute viscosity (uo) 0.138 Ns/m? 
Modulus of elasticity (E,, E2) 20-72 1082 N/m? 
Borsson’ si ratio. (v7) 05) 0.3 

Cylinder rotating speeds (u;, uz) 3.048 m/s 
Pressure viscosity exponent (a) ts 2 hoe a/N 
Minimum film thickness (h,,;,) Pees xtOr Asam to 


SUG xt Ome. mm 
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